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Using an effective low-energy Hamiltonian derived from the first-principles 
electronic structure calculations for the narrow t2 g bands of YTiC>3, LaTi03, 
YVO3, and LaVC>3, we evaluate the contributions of the correlation energy 
(Ec) to the stability of different magnetic structures, which can be realized 
in these distorted perovskite oxides. We consider two approximations for Ec, 
which are based on the regular perturbation theory expansion around a nonde- 
generate Hartree-Fock ground state. One is the second order of perturbation 
theory, which allows us to compare the effects of local and nonlocal correla- 
tions. Another one is the local t-matrix approach, which allows us to treat 
some higher-order contributions to Ec- The correlation effects systematically 
improve the agreement with the experimental data and additionally stabilize 
the experimentally observed G- and C-type antiferromagnetic (AFM) struc- 
tures in YVO3 and LaV03, though the absolute magnitude of the stabilization 
energy is sensitive to the level of approximations and somewhat smaller in the 
t-matrix method. The nonlocal correlations additionally stabilize the ferro- 
magnetic ground state in YTiC>3 and the C-type AFM ground state in LaV03. 
Amongst two inequivalent transition-metal sites in the monoclinic structure, 
the local correlations are stronger at the sites with the least distorted environ- 
ment. Limitations of the regular perturbation-theory expansion for LaTiC>3 
are also discussed. 

PACS: 71.10.-w; 71.15.Nc; 71.28.+d; 75.25.+Z 
1. Introduction 

An interest to the transition-metal perovskite oxides YTiC>3, LaTiC>3, YVO3, and 
LaVC>3 is mainly related with the variety of magnetic structures, which can be real- 
ized in these, seemingly alike, compounds. For example, YTiC>3 has the ferromag- 
netic (F) structure pQ. LaTiC>3 is a three-dimensional (G-type) antiferromagnet [2]. 
At the low temperature, YVO3 forms the G-type antiferromagnetic (AFM) struc- 
ture, which can be transformed to a chainlike (G-type) antiferromagnetic structure 
at around 77 K [3]. On the contrary, LaV03 is the G-type antiferromagnet in the 
whole temperature range below the magnetic transition temperature j3j. Surpris- 
ingly, the difference exists not only between titanates (YTi03 and LaTiOs) and 
vanadites (YVO3 and LaVOs), which have a different number of valent electrons, 
but also within each group of formally isoelectronic materials. The differences are 
apparently related with the tiny changes in the distorted perovskite structure, which 
are amplified by the effects of Coulomb correlations in the narrow t2 9 band. The 
details of the crystal structure can be found in Refs. [TJ 12 EJ H|- Briefly, both ti- 
tanites have an orthorhombic structure, although the details of this structure are 
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rather different for YTiOa and LaTi03. LaV03 is crystallized in a monoclinic struc- 
ture. The low-temperature phase of YVO3 is orthorhombic (shown in Fig. QJ, which 
becomes monoclinic at around 77 K. The structural orthorhombic-monoclinic tran- 
sition coincides with the G-C AFM transition. Generally, Y-based compounds are 
more distorted (due to smaller size of the Y 3+ ions). 

There is a large number of theoretical articles devoted to the origin of the mag- 
netic ground states in the distorted t 2g perovskite oxides. The problem has been con- 
sidered on the basis of first-principles electronic structure calculations (e.g., Refs. 
and the model approaches for the strongly-correlated systems (e.g., Refs. [E10IE])- 
The model theories typically vary on the assessment of the role played by the lattice 
distortions and the Coulomb correlations |H] . 

We believe that any realistic theoretical description of these compounds is prac- 
tically impossible without the impact from the first-principles electronic structure 
calculations: simply, the lattice distortion is too complex, and, had we try to pos- 
tulate a model Hamiltonian for these t 2g perovskite oxides, we would inevitably 
face the problem of choosing the values for a large number of model parameters, 
which cannot be fixed in unbiased way. However, the conventional electronic struc- 
ture calculations are also far from being perfect. Typically, they are supplemented 
with some additional approximations, which have serious limitations for treating the 
Coulomb correlations in the case of strongly-correlated materials. A typical exam- 
ple is the local-density approximation (LDA). From this point of view, a promising 
direction is to make a bridge between first-principles electronic structure calcula- 
tions and models for the strongly-correlated systems, and construct an appropriate 
model Hamiltonian entirely "from the first principles" . Fortunately, in the case of 
transition-metal oxides, we are typically dealing only with a small group of states 
located near the Fermi level and well separated from the remaining part of the spec- 
trum (for instance, t 2g bands in Fig. [TJ. These states are mainly responsible for 
the electronic and magnetic properties of oxide materials. Therefore, in may cases 
it is sufficient to consider a minimal model, consisting of only the t 2g bands, and 
include the effect of other bands into the renormalization of interaction parameters 
in the t 2g band. Such a strategy was pursued in Refs. [HI EH HI] ■ It consists of three 
major steps: first-principles electronic structure calculations — > construction of the 
model Hamiltonian — > solution of this model Hamiltonian. The first applications 
to the distorted t 2g perovskite oxides have been considered in Refs. [TTJ H2] . The 
present paper deals with the last part of the problem. We will solve the model 
Hamiltonian derived in Ref. jllj . and mainly focus on the role played by the correla- 
tion effects, beyond the mean-field Hartree-Fock (HF) approximation. Particularly, 
we will consider two perturbative approaches. One is the regular second-order per- 
turbation theory for the correlation energy [T^], and the other one is the t-matrix 
approach jT3J [TSJ HH] . In both approaches, the HF approximation is used as the 
starting point. This implies that the degeneracy of the HF ground state is already 
lifted by the crystal distortion so that the regular perturbation theory is justified. 
We will also discuss some limitations of this treatment for LaTi03. 
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2. Construction of model Hamiltonian 

Our first goal is the construction of the effective multi-orbital Hubbard model for 
the isolated t 2g bands: 

^ = X! X! + 9 X! X! ^a/37<^Ra^R7^R/^RS' ( X ) 

RR' a/3 R a/3j5 

where c Ra (crq,) creates (annihilates) an electron in the Wannier orbital of 
the transition-metal site R, and a is a joint index, incorporating all remaining (spin 
and orbital) degrees of freedom. The matrix ^rr' = |I^rr.'II parameterizes the kinetic 
energy of electrons, where the site-diagonal part (R=R/) describes the local level- 
splitting, caused by the crystal field, and the off-diagonal part (R^R') stands for 
the transfer integrals. U aPjS = J dr J dr'W£\r)Wl(r)v SCI (r-r')Wg (r')Wi(r') = 
(a7|v scr |/55) are the matrix elements of screened Coulomb interaction v scr (r— r'), 
which are supposed to be diagonal with respect to the site indices. In principles, 
can also depend on the site-index R. Nevertheless, for the sake of simplicity 
of our notations, here and throughout in this paper we drop the index R in the 
notation of the Coulomb matrix elements (however, we do consider this dependence 
in all our calculations). 

The procedure of mapping of the first-principles electronic structure calculations 
onto the model Hamiltonian (0) for distorted perovskite oxides has been discussed 
in details in Refs. [TO! lllj. Here, we only outline the main idea. The kinetic-energy 
part, /irr', can be obtained using the downfolding method, which is exact and 
equivalent to the projector-operator method ^7j. The Wannier functions can be 
formally derived from /irr/, using the ideas of the linear-muffin-tin-orbital (LMTO) 
method fOJ The matrix of screened Coulomb interactions in the t2 9 band 

can be calculated using a hybrid approach, which combines the constraint density- 
functional theory with the random-phase approximation for the hybridization effects 
between transition-metal d and other atomic states |1U| . The values of the model 
parameters obtained in such a way can be found in Ref. [TT] . 

3. Solution of model Hamiltonian 

3.1. Hartree-Fock approximation 

The HF method provides the simplest approximation for the solution of the many- 
electron problem with the Hamiltonian (0). In this case, the trial many-electron 
wavefunction is searched in the form of a single Slater determinant \S{(p k }), con- 
structed from the one-electron orbitals {y?£,}. In these notations, k is a joint index, 
which contains the information about the momentum (k) in the first Brillouin zone, 
the number of band (n), and the spin (a= f or j) of the particle. The one-electron 
orbitals {tpk} are subjected to the variational principle and requested to minimize 
the total energy 

E HF = min(S{(p k }\H\S{(p k }) 
for a given number of particles Af. This yields the following equations for {f k }- 

(h k + $)\(p k )=e k \<p k ), (2) 
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where /ik=||^k^|| is the kinetic part of the model Hamiltonian in the reciprocal 
space: h^f=j] J2w ^RR' e_ik '( R ~ R ') (AT being the number of sites), and V"=||V^^|| is 
the HF potential: 

Va/3 — 12 (Ua/3-yS ~ U a s~,p) Tl^g- (3) 

In the following, we will also use the notation /ihf> which stands for the total Hamil- 
tonian of the HF method, h+V. Eq. (|2j) is solved self-consistently together with the 
equation 

occ 

" = 12 lv 3 fc)(v ? fc| 

k 

for the density matrix n^Hn^H in the basis of Wannier orbitals. Finally, the total 
energy in the HF method can be obtained as 

OCC ^ 

-^HF = 12 £k ~ o 12 V/3oin>aP- 
k Z a/3 

3.2. Second Order Perturbation Theory for Correlation Energy 

The simplest way to go beyond the HF approximation is to include the correlation 
interactions in the second order of perturbation theory for the total energy [T^]. The 
correlation interaction (or a fluctuation) is defined as the difference between true 
many-body Hamiltonian and its one-electron counterpart, obtained at the level 
of HF approximation: 

= 12 f \ 12 ^a/3 7 5cLcR 7 C R/3 C R(5 - £ V^C^J^ j . (4) 
R \ Z af3y8 a/3 j 

By treating He as a perturbation, the correlation energy can be easily estimated 
as [El: 

E (2) _ ( G \7~Lc\S) (S\H C \G) ^ 

where \G) and IS*) are the Slater determinants corresponding to the low-energy 
ground state in the HF approximation, and the excited state, respectively. Due to 
the variational properties of the HF method, the only processes which may con- 
tribute to E c are the two-particle excitations, for which each IS") is obtained from 
\G) by replacing two one-electron orbitals, say tp^ and tp k ^ from the occupied part 
of the spectrum by two unoccupied orbitals, say ^ fc3 and tpk 4 - Hence, using the 
notations of Sec. the matrix elements take the following form: 

(S\H C \G) = (hk^v^k^) - (k 3 k4\v SCI \k2h) . (6) 

Then, we employ a common approximation of noninteracting quasiparticles and 
replace the denominator of Eq. (jSJ) by the linear combination of HF eigenvalues: 
Euf(S)— E-ft-p(G) ~ Sk z +Sk i —£k 1 —^k2 [XHJ - The matrix elements (jHJ) satisfy the 
following condition: (S\Hc\G)~jj J2n e i(k3+k4 ~ kl ~ k2) ' R , provided that the screened 
Coulomb interactions are diagonal with respect to the site indices. 

A good point of the second-order of the perturbation theory is that it allows 
us to estimate relatively easily both on-site (R=0) and intersite (R^O) elements 
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of this expansion. In the following, we will use this method in order to study the 
relative role played by these effects in the stability of different magnetic structures 
of the distorted perovskite oxides. The R=0 term corresponds to the commonly 
used single-site approximation for the correlation interactions, which becomes exact 
in the limit of infinite spacial dimensions [T9*] . 



3.3. t-matrix approach 

The basic idea of the t-matrix approach is to look at the true many-electron system 
as a superposition of independent two-electron subsystems, and to solve rigorously 
the Schrodinger equations for each of these subsystems [HI UHl EE]. Hence, we 
consider the following two-electron Hamiltonian: 

H(l, 2) = ^hf(I) + h BF (2) + Av(l, 2), 

where Av(l, 2)—v scr (l, 2)— V^l)— V(2), v scr (l,2) is the screened (by other bands) 
Coulomb interactions between electrons '1' and '2' in the t 2g band, and h^F {V) is 
the one-electron Hamiltonian (potential) in the HF approximation. For the periodic 
system, the Schrodinger equation can be written in the following form: 

H\^f kl k 2 ) = E klk2 \ty klk2 }. (7) 

Any two-electron wavefunction \^ kl k 2 ) can be expanded in the basis of (also two- 
electron) Slater's determinants: \kik 2 ) = -^{(p kl (l)(p k2 (2) — ^ (1)^(2)}, etc. Apart 
from a normalization multiplier, this expansion has the following form |16j : 

|** lfa > = M»>+ E r£&IM*>- (8) 
\k3k4) 

Note that the summation goes only over nonequivalent Slater's determinant 
\k3k4), constructed from the one-electron orbitals k% and k&. For example, since 
I ^4^3) = — \k3k4), the determinant \k4k3) should be excluded from the sum (jHJ), 
etc. By substituting Eq. (jSJ) into Eq. (JZJ), and introducing the new notations 
^E klk2 =E klk2 —e kl —e k2} such that 

h UF (l) + h n¥ (2) - e kl - e fc J \k x k 2 ) = 

(i.e., e kl and e k2 are the eigenvalues of the HF Hamiltonian), one obtains the follow- 
ing equation for AE klk2 and r fc \*: 

(Av - AE klk2 ) \k Y k 2 ) + J2 ( £ k 3 + ^ - £fc! - ?k 2 +Av- AE klk2 ) T k k f k \\hk A ) = 0. 

\k s k 4 ) 

By considering the matrix element of this equation with {kik 2 \, one can find that 
AE klk2 = (hk^Avlhk,) + J2 T^ihhlAvlhh), (9) 

where the first term is the energy of Coulomb and exchange interactions in the HF 
approximation (minus the potential energy), while the second term is the correlation 
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energy. By considering the matrix elements with (A^Afel^fciA^I, one can find another 
set of equations for r^ fc *: 

(A: 5 A; 6 |A^|A; 1 fc 2 ) + (e A;5 + - s kl - e k2 - AE klk2 ) rjg£+ £ ^(^elA^M = 0. 

\k3k4) 

They are solved iteratively with respect to Av. In order to do so, it is convenient 
to introduce the two-particle Green's function, 

g,^ k _ \hh){hh\ 

\k z ki) £fc 3 + £ ki - - ^fc 2 _ AE klk2 ' 

and derive a matrix equation for {T^*}, which are then substituted into Eq. Q. 
Then, it is rather straightforward to derive the following expression for AE klk2 : 

AE klk2 = {k 1 k 2 \T kxk2 \kik 2 ) , (10) 

where T klk2 is the so-called t-matrix: 



T klk2 = Av 



l + G klk2 Av . (11) 



The correlation energy of the t-matrix method is obtained after the subtraction from 
Eq. (|10|) the energies of Coulomb and exchange interactions in the HF approximation 
and summation up over all Slater's determinants constructed from the occupied one- 
electron orbitals of the HF method: 

occ 

E c = E (kik 2 \f klk2 - Avlhh). (12) 

\kik 2 ) 

In practice, each HF orbital has been expanded over the basis of Wannier functions, 
and then all calculations of T klk2 and Eq have been performed in this basis. 

By expanding T klk2 up to the second order of Av, we regain Eq. (JHJ), obtained 
in the second order of perturbation theory. Therefore, the good point of the t- 
matrix approach is that it allows us to go beyond the second order of perturbation 
theory and evaluate the higher order effects of Av onto the correlations energy. 
Nevertheless, it was supplemented with some additional approximations. 

1. When we compute the matrix elements of the form (^3^4! Av\kik2), being 
proportional to -h J2n e^ k3+k4 ~ kl ~ k2 ^' R , we consider only the R=0 part of this 
sum and neglect all other contributions. This corresponds to the single-site 
approximation for the t-matrix. 

2. In all matrix elements (k^k^ Av\kik2), we replace Av by v SCT and drop the 
one-electron potentials of the HF method. Strictly speaking, this procedure 
is justified only when both one-electron states k\ and k 2 are different from 
k 3 and fc 4 , for example, when they belong, correspondingly, to the occupied 
and unoccupied part of the spectrum, like in the second order of perturbation 
theory. However, this is no longer true for the higher-order terms with respect 
to Av. Nevertheless, we believe that the difference is small. 

All correlation energies have been computed in the mesh of 75 points in the first 
Brillouin zone (BZ), corresponding to the 4:4:2 divisions of the reciprocal translation 
vectors for the distorted perovskite structure. The actual integration over the BZ 
has been replaced by the summation over this mesh of points. 
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4. Results and Discussions 

First applications of the proposed method to YTiC>3, LaTiOs, YVO3, and LaVC>3 
have been considered in Ref. [TT], where we have summarized results of HF cal- 
culations for the model (0) and the behavior of correlation energies in the second 
order of perturbation theory, supplemented with the single-site approximation. In 
the present work we will further elaborate the problem by focusing on the following 
question: 

1. the role of higher-order contributions to the correlation energy; 

2. the role of nonlocal (or intersite) contributions to the correlation energy 

We will also consider the effects of monoclinic distortion and analyze the contribu- 
tions to the correlation energy of inequivalent transition-metal sites. The results of 
these calculations are presented in Tables HE for all considered compounds. First, 
we would like to summarize the main results of Ref. JT] . 

1. The HF approximation yields the correct magnetic ground state for YTiC>3, 
LaV03, and both phases of YVO3. This conclusion is fully consistent with 
the results of accurate all-electron band-structure calculations [Sj, and this is 
quite remarkable that all these results can be reproduced in our minimal model 
derived for the t 2g bands. 

2. The correlation effects favor the AFM spin alignment and additionally stabilize 
the experimentally observed G- and C-type AFM states in YVO3 and LaV03. 

3. None of the considered approaches reproduces the experimental G-type AFM 
ground state of LaTiOs (instead, the theoretical calculations steadily converge 
to the A-type AFM ground state [HUE!)- 

Then, what will happen if we go beyond the second-order perturbation theory and 
apply the t-matrix approach? Generally, the t-matrix approach reduces the absolute 
value of the correlation energy. However, the magnitude of this reduction strongly 
depends on the magnetic state. For example, if the F state is only weakly affected 
by the higher-order correlation effects (the typical changes of Eq varies from 1% in 
YVO3 till 13% in LaTiO-3), E c in the G-type AFM phase can drop by nearly 50%. 
From this point of view, if the second order or perturbation theory does not solve the 
problem of the G-type AFM ground state of LaTi03, it seems to be unlikely that the 
higher-order effects can reverse the situation. Apparently, LaTi03 is different from 
other perovskite oxides, and the regular perturbation-theory expansion, though may 
be justified for the majority of considered compounds, does not work in the case of 
LaTiO"3. This seems to be reasonable, because LaTiO"3 has the largest correlation 
energies, which are comparable with the splitting of the t<i g levels caused by the 
crystal distortion (~37 meV |Hj). Therefore, it is quite possible that the correlation 
effects in LaTiO"3 should be considered at the first place, and the simple HF theory 
for the spin and orbital ordering with the subsequent inclusion of the correlation 
effects as a perturbation to the HF ground state may not be appropriate here IB] . 
Note that in other materials, the situation is different: the typical values of the t2 9 - 
levels splitting in YTiC>3, YVO3, and LaVOs are about 100 meV [11J, which exceeds 
the correlation energy by at least one order of magnitude. Therefore, it seems that 
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the degeneracy of the HF ground state is already lifted by the crystal distortion, 
and the correlation effects are well described by means of the regular perturbation 
theory expansion. This is partly supported by recent total energy calculations for 
the orthorhombic phase of YVO3 using path-integral renormalization group method, 
which is free of any perturbation-theory expansions for the correlation energy [20] • 
The method was applied to the same model, and the main conclusions concerning 
the magnetic phase diagram were similar to our present finding. 

The correlations additionally stabilize the experimentally observed C- and C- 
type AFM states in YVO3 and LaVOs. Moreover, in the orthorhombic phase of 
YVO3, the correlation effects tend to stabilize the C-type AFM state; while in the 
monoclinic phase, they stabilize the C-type AFM state, being in total agreement 
with the experimental data. This trend is clearly seen both in the second order of 
perturbation theory and in the t-matrix approach, though the latter yield somewhat 
smaller values for the stabilization energy associated with the correlation effects. 

The higher-order correlations play an important role YTiC>3 and additionally 
stabilize the ferromagnetic phase. The latter emerges as the ground state already in 
the HF approach, where the total energy difference between ferromagnetic and the 
next A-type AFM state is about 2.05 meV per one Ti site (Table GJ. However, if we 
take into account the correlation effects in the second order of perturbation theory 
(and consider the single-site approximation), this difference is reduced to only 0.99 
meV. Therefore, the situation is very fragile. Nevertheless, the t-matrix approach, 
which affects more strongly the A-type AFM state, will recover some of these energy 
gains and make the total energy difference between ferromagnetic and A-type AFM 
states to be about 1.83 meV per one Ti site. 

The intersite correlation energies, which have been estimated in the second order 
of perturbation theory, can be large in some ferromagnetically couple bonds. This 
is especially true for YTiC>3 and LaVC>3. For example, the energy of interaction 
between nearest-neighbor sites '1' and '2' (see Fig.QJ, located in the ab-plane of the 
ferromagnetic phase of YTiC>3 is about —0.38 meV (TableCJ). Since in the ab-plane, 
each transition-metal atom interacts with four nearest neighbors, it corresponds to 
the additional energy gain —0.38x4 = —1.52 meV per one Ti site. Similar estimates 
yield -0.37x4 = -1.48 meV, -0.17x4 = -0.68 meV, and -0.18x4 = -0.72 meV, 
correspondingly for the A-, C-, and C-type AFM states. Therefore, the in-plane 
intersite correlations tend to additionally stabilize the ferromagnetic phase relative 
to the AFM states C and G. In the A-type AFM phase, the sites 'l'and '2' are also 
ferromagnetically coupled, like in the totally ferromagnetic phase. Therefore, these 
two phases have practically the same intersite correlation energies in the ab-plane. 
The inter-plane correlations appear to be small in all magnetic phases of YTiCV 

In LaVC>3, the situation is somewhat different, and this is a good example of 
the system where already the inter-plane correlations play a more important role. 
Indeed, the energies of intersite correlations are the largest in the ferromagnetic 
chains of C-type AFM phase, which is also the magnetic ground state of this com- 
pound. These energies are associated with the bonds '1-3' and '2-4', which are 
shown in Fig. Q and the results are summarized in Table EJ Thus, in the case of 
LaV03, the inter-plane correlations additionally stabilize the C-type AFM ground 
state. However, since each transition-metal atom interacts with only two nearest 
neighbors along the c axis, the stabilization energy is not particularly large (about 
—0.26x2 = —0.52 meV per one V atom). 
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The intersite correlation energies are large also in the case of LaTi03 (Table EJ). 
However, they tend to stabilize either ferromagnetic or A-type AFM states, and do 
not explain the appearance of experimental G-type AFM ground state. Again, we 
believe that the problem is related with the use of the regular perturbation theory 
expansion, which may not be justified in the case of LaTiCV 

The monoclinic distortion realized in LaVOs and in the high-temperature phase 
of YVO3 produces two inequivalent pairs of transition-metal sites, which are shown 
correspondingly as (1,2) and (3,4) in Fig. ^ Therefore, it is interesting to consider 
the interplay between correlation energies and the lattice distortions around differ- 
ent transition-metal sites. In our notations, the crystal structure around the sites '3' 
and '4' is more distorted than the one around the sites '1' and '2'. Such a distortion 
directly correlates with the magnitude of the crystal-field splitting in different sub- 
lattices jllj . Then, the on-site correlations are generally stronger at the sites with 
the least distorted environment (site '1' in the Tables El and EJ). This rule holds both 
for YVO3 and LaV03 (though with some exception for the ferromagnetic phase of 
LaVOs). I n the C-type AFM phase, which is always realized as the magnetic ground 
state in the monoclinic structure, the difference of on-site correlation energies asso- 
ciated with the sites '1' and '4' is about 1 meV per one V site, as obtained in the 
second order of perturbation theory. This value is further reduced till 0.5 meV per 
one V site by higher-order correlations in the t-matrix theory. 

5. Summary and Conclusions 

This paper is the continuation of previous works (Refs. [Ej, [H], and [12]) devoted 
to the construction and solution of an effective low-energy models for the series of 
distorted ti g perovskite oxides on the basis of first-principles electronic structure 
calculations. It deals with the analysis of correlation interactions and their contri- 
butions to stability on different magnetic structures, which can be realized in these 
compounds. The correlation energies have been calculated on the basis of a regular 
perturbation theory expansion starting from the ground state of the HF method. 
Thus, our strategy implies that the degeneracy of the HF ground state is already 
lifted by the crystal distortion and the regular perturbation theory is justified. This 
seems to be a good approximation for the most distorted YTiC>3, YVO3, and even 
LaVC>3, where 

1. the correct magnetic ground state can be formally obtained at the level of HF 
approximation; 

2. the correlation effects, included as a perturbation to the HF ground state, 
systematically improve the agreement with the experimental data. 

However, in LaTi03, the situation is completely different: 

1. the HF method yields an incorrect magnetic ground state (A- type AFM in- 
stead of G-type AFM); 

2. the correlation interactions, treated as a perturbation to this incorrect HF 
ground state, do not change the overall picture, and the G-type AFM state 
remains unstable relative to the A state. 
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Thus, the origin of the G-type AFM ground state in LaTiOs seems to be different 
from other perovskite oxides and remains a challenging problem for future theories. 
Apparently, one of our basic assumptions about the nondegeneracy of the HF ground 
state breaks down in the case of LaTi03, and the true ground state cannot be 
approached through the series of continuous corrections applied to the single-Slater- 
determinant HF theory. Therefore, the next important step for LaTiOs would be 
get rid of this "nondegeneracy assumption" and expand the class of the possible 
ground states, which would include some aspects of the orbital liquid theory [Jj. 
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Table 1: Hartree-Fock, -Ehf, and correlation energies obtained in the second order of 
perturbation theory, Eq\ and in the t-matrix approach, Eq \ for the orthorhombic 
phase of YTiC>3. The Hartree-Fock energies are measured from the most stable 
magnetic state in meV per one formula unit. The correlation energies are measured 
in meV per one transition-metal site or a pair of sites, correspondingly for the on-site 
and intersite contributions. Note that the t-matrix was computed in the single-site 
approximation. Therefore, only the site-diagonal part of Eq* is shown. The positions 
of the transition-metal sites are shown in Fig. ^ 



Phase 


Ebf 








E (t) 






Ti! 


Ti r Ti 2 


Tii-Ti 3 


Tii 


F 





-5.13 


-0.38 


-0.01 


-4.58 


A 


2.05 


-6.19 


-0.37 





-4.80 


C 


14.40 


-8.32 


-0.17 


-0.01 


-5.28 


G 


16.25 


-8.48 


-0.18 


-0.01 


-5.31 
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Table 2: Hartree-Fock, Ebf, and correlation energies obtained in the second or- 
der of perturbation theory, Eq \ and in the i-matrix approach, Eq , for the low- 
temperature orthorhombic phase of YVO3 (T< 77 K). All energies are measured in 
meV. See Fig. ^for the details of the notations. 



Phase 


Euf 




4 2) 




Eg) 






Vi 


Vi-V 2 


v r v 3 


Vi 


F 


21.66 


-2.19 


-0.12 


-0.02 


-2.16 


A 


14.59 


-4.67 


-0.12 


-0.01 


-3.31 


C 


10.14 


-5.61 


-0.07 





-3.14 


G 





-7.07 


-0.07 


-0.01 


-4.06 
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Table 3: Hartree-Fock, Ebf, and correlation energies obtained in the second or- 
der of perturbation theory, Eq\ and in the t- matrix approach, Eq , for the high- 
temperature monoclinic phase of YVO3 (77 K <T< 116 K). All energies are mea- 
sured in meV. See Tabled for the details of the notations. Note that, in the mon- 
oclinic phase, the planes 1-2 and 3-4 are inequivalent (see Fig. Q). Therefore, there 
are two different types of on-site (denoted as Vi and V 4 ) and intersite (denoted as 
V1-V2 and V4-V3) contributions to the correlation energy. The contributions V1-V3 
and V4-V2 are equivalent and are both shown only for the sake of completeness. 



Phase 


Erf 






4 2) 








E {t) 








Vi 


V1-V2 


V1-V3 


v 4 


V4-V3 


V4-V2 


Vx 


v 4 


F 


11.71 


-2.81 


-0.02 


-0.03 - 


-1.74 


-0.01 


-0.03 


-2.76 - 


-1.71 


A 


13.97 


-5.87 


-0.03 


-0.01 - 


-3.63 


-0.01 


-0.01 


-4.14 - 


-2.55 


C 





-8.08 


-0.02 


-0.05 - 


-6.98 


-0.03 


-0.05 


-4.85 - 


-4.33 


G 


6.63 


-7.56 


-0.02 


-0.01 - 


-6.49 


-0.03 


-0.01 


-4.38 - 


-3.76 
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Table 4: Hartree-Fock, Ehf, and correlation energies obtained in the second order 
of perturbation theory, Eq\ and in the t- matrix approach, Eq , for the monoclinic 
phase of LaVOs. All energies are measured in meV. See Tables ^ and H3 for the 
details of the notations. Note that, in the monoclinic phase, the planes 1-2 and 
3-4 (see Fig. [TJ are inequivalent, that results in two types of V sites as well as the 
in-plane interactions. 



Phase 


Erf 






















Vi 


v r v 2 


Vi-V 3 


v 4 


V4-V3 


v 4 -v 2 




v 4 


F 


20.98 


- 3.82 


-0.02 


-0.15 - 


4.13 


-0.02 


-0.15 


-3.74 - 


-4.02 


A 


20.63 


-11.77 


-0.22 


-0.03 - 


8.80 


-0.02 


-0.03 


-8.34 - 


-5.84 


C 





-13.37 


-0.04 


-0.26 - 


12.54 


-0.02 


-0.26 


-8.86 - 


-8.39 


G 


7.65 


-10.52 


-0.04 


-0.02 - 


9.02 


-0.03 


-0.02 


-6.17 - 


-5.41 
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Table 5: Hartree-Fock, £7hf, and correlation energies obtained in the second order of 
perturbation theory, Eq , and in the t-matrix approach, Ej& , for the orthorhombic 
phase of LaTiOg. All energies are measured in meV. See Tabled for the details of 
the notations. 



Phase 










h c 






Tlx 


Ti r Ti 2 


Ti r Ti 3 


Tii 


F 


4.95 


-11.08 


-0.52 


-0.08 


- 9.66 


A 





-22.53 


-0.54 


-0.07 


-15.17 


C 


19.57 


-17.19 


-0.23 


-0.11 


-11.04 


G 


11.51 


-23.02 


-0.22 


-0.09 


-13.99 
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Figure 1: A characteristic example of the crystal structure (left) and the electronic struc- 
ture in the local-density approximation (right) of the orthorhombically distorted YVO3. 
In the left panel, the symbols a, b, and c stand for orthorhombic translations, and the 
symbols 1-4 denote the transition- metal sites, which form the unit cell of the distorted 
perovskite oxides. In the right panel, the shaded area shows contributions of the atomic 
V(3d) states. Other symbols show the positions of the main bands. The Fermi level is at 
zero energy. 
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